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ABSTRACT 

We report our measurement of the rate of change of period with time (P) for 
the 215 s periodicity in the pulsating white dwarf G 117-B15A, the most stable 
optical clock known. After 31 years of observations, we have finally obtained a 
4(7 measurement Pobserved = (4.27 ± 0.80) x lO^^^s/s. Taking into account the 
proper-motion effect of Pproper = (7.0±2.0) x 10^^^ s/s, we obtain a rate of change 
of period with time of P = (3.57 ± 0.82) x 10~^^s/s. This value is consistent 
with the cooling rate in our white dwarf models only for cores of C or C/0. With 
the refinement of the models, the observed rate of period change can be used to 
accurately measure the ratio of C/0 in the core of the white dwarf. 

Subject headings: Stars: evolution - stars: oscillations - stars: individual: G 117- 
B15A 

1. Introduction 

G 117-B15A is a pulsating white dwarf with a hydrogen atmosphere, a DAV, also called 
a ZZ Ceti star (McGraw 1979). These stars show multi-periodic non-radial g-mode pulsations 
that can be used to measure their internal properties and rate of evolution. 
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McGraw & Robinson (1976) found the star was variable, and Kepler et al. (1982) studied 
its light curve, finding 6 simultaneous pulsations. The dominant mode has a period of 215 s, 
a fractional amplitude of 22 mma (milli- modulation amplitude = 1/1.086 milli magnitude), 
and is stable in amplitude and phase. The other smaller pulsation modes vary in amplitude 
from night to night (Kepler et al. 1995), suggesting the presence of unresolved components. 
Because the DAVs appear to be normal stars except for their variability (Robinson 1979; 
Bergeron et al. 1995, 2004), i.e., an evolutionary stage in the coohng of all white dwarfs, it 
is likely that the DAV structural properties are representative of all DA white dwarfs. 

Mukadam ct al. (2004a) and MuUally et al. (2005) discovered 46 new ZZ Cctis and 
more than doubled the number of known variables, using Tcs and log g values derived from 
the optical spectra obtained by the Sloan Digital Sky Survey (Kleinman et al. 2004). While 
Bergeron et al. (1995, 2004) find the ZZ Ceti instability strip to be pure, i.e., contain no non- 
variable stars, Mukadam et al. (2004a,b); MuUally et al. (2005) find several stars inside the 
instability strip for which they could detect no variability above their detection threshold. It 
is necessary to obtain S/N > 50 spectra of the non- variables within the strip and rc-analyze 
their Tes and log g values, and also to obtain additional time series photometry on these 
stars, to ensure that they are not low amplitude variables, and re-study the purity of the 
instability strip. 

We report our continuing study of the star G 117-B15A, also called RY LMi and 
WD 0921+354, one of the hottest of the ZZ Ceti stars. We expect that the rate of change 
of a pulsation period with time for (yf-modc pulsations in white dwarf stars to be directly 
related to its evolutionary timescale (Winget et al. 1983), allowing us to infer the age of a 
cool white dwarf since its formation. We have been observing the star since 1974 to measure 
the rate of period change with time (P) for the largest amplitude periodicity. Using all the 
data obtained from 1974 through 1999, Kepler et al. (2000) arrived at a determination of 
P= (2.3 ±1.4) X 10-^5 s/s. 

G 117-B15A was the first pulsating white dwarf to have its main pulsation mode index 
identified. The 215 s mode is an ^ = 1, as determined by comparing the ultraviolet pulsation 
amphtude (measured with the Hubble Space Telescope) to the optical amphtude (Robinson 
et al. 1995). Kotak et al. (2004) confirm the i measurement for the P=215 s pulsation, and 
show that the other large amplitude modes, at 271 s and 304 s, show chromatic amplitude 
changes that do not fit the theoretical models, using time-resolved spectra obtained at the 
Keck Telescope. Robinson et al. (1995), and Kocstcr, Allard, & Vauclair (1994) derive T^s 
near 12,400 K, while Bergeron et al. (1995, 2004) using a less efficient model for convection, 
derive Tefr=ll,630 K. 

Kepler (1984) demonstrated that the observed variations in the light curve of G 117- 
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B15A are due to non-radial ^mode pulsations. Kepler et al. (2000) show the models predict 
the effect of radius change due to still ongoing contraction to be an order of magnitude 
smaller than the cooling effect on the rate of period change. 

G 117-B15A is proving to be a useful laboratory for particle physics (Isern et al. 2004). 
Corsico et al. (2001) calculated the hmit on the axion mass compatible with the then observed 
upper limit to the cooling, showing maCOS/3 < 4.4 meV and Kepler (2004) demonstrates 
axion cooling would be dominant over neutrino cooling for the lukewarm white dwarf stars 
for axion masses of this order. Biesiada & Malec (2002) show that the 2(7 upper limit 
published in Kepler et al. (2000) limits the string mass scale Af^ > 14.3 TeV/c^ for 6 
dimensions, from the observed cooling rate and the emission of Kaluza-Klein gravitons, but 
the limit is negligible for higher dimensions. Benvenuto et al. (2004) show the observed rates 
of period change can also be used to constrain the dynamical rate of change of the constant 
of gravity G. 

Bradley (1996, 1998) used the mode identification and the observed periods of the 3 
largest known pulsation modes to derive a hydrogen layer mass lower limit of 10~^M*, and 
a best estimate of 1.5 x 10~^M*, assuming A; = 2 for the 215 s mode, and 20:80 C/0 core 
mass. The core composition is constrained mainly by the presence of the 304 s pulsation. 
Benvenuto et al. (2002) show the seismological models with time-dependent element diffusion 
are only consistent with the spectroscopic data if the modes are the 1=1^ k=2, 3, and 4, 
and deduces M = 0.525 Mq, logiMn/M^) > -3.83 and T^s = 11800 K. Their best model 
predicts: parallax n=15.89 mas, P = 4.43 x 10"^^ s/s, for the P=215 s, P = 3.22 x 10"^^ s/s, 
for the P=271 s, and P = 5.76 x 10~^^ s/s, for the P=304 s periodicities. 

2. Observations 

Kepler et al. (2000) reported on the observations from 1974 to 2000. We report in this 
paper additional 19.3 h of time series photometry in 2001, 30.6 h in 2002, 24.4 h in 2003, 
4 h in 2004, and 13.6 h in 2005, most using the Argos prime-focus CCD camera (Nather & 
Mukadam 2004) on the 2.1 m Otto Struve telescope at McDonald Observatory. 

We observed the light through a BG40 filter to maximize the signal-to-noise ratio (S/N) 
because the pulsation amplitudes are small (2%), the star is faint (V=15.52, Eggen & Green- 
stein 1965), and also because the non-radial g^mode light variations have the same phase in 
all colors (Robinson, Kepler, & Nather 1982) but the amplitudes decrease with wavelength. 
For example, a filter-less observation with Argos gives an amplitude around 40% smaller for 
G 117-B15A. 
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3. Data Reduction 

We reduce and analyze the data in the manner described by Nather et al. (1990), and 
Kepler (1993). We bring all the data to the same fractional amplitude scale, and the times 
from UTC to the uniform Barycentric Julian Coordinated Date (TCB) scale, using JPL 
DE405 ephemeris (Standish 1998, 2004) to model Earth's motion. We compute Fourier 
transforms for each individual run, and verify that the main pulsation at 215 s dominates 
each data set and has an amplitude stable up to 15%, our uncertainty in amplitude due to 
the lack of accurate time- and color-dependent extinction determination. 



4. Time Scale for Period Change 

As the dominant pulsation mode at P=215 s has a stable frequency and amplitude since 
our first observations in 1974, we can calculate the time of maximum for each new run and 
look for deviations. 

We fit our observed time of maximum light to the equation: 

(O-C)^ AEo + AP-E+-P-P-E^ 

where AEo = {T^ax - ^max)> AP = (P - Pt=T^^J, E is the epoch of the time of maximum, 
i.e, the integer number of cycles after our first observation, T^gj, is the time of maximum 
assumed, T^^^^ is the time of maximum that best fit the parabola, P is the period that best 
fit the parabola, and Pt=T°^^^ is the period assumed at T^^^ ^. The times of maxima are 
calculated by a linear least-squares fit of the light curve of each night to a sum of the six 
detected frequencies. They are shown in Table 6. 

In Figure 1, we show the 0-C timings after subtracting the correction to period and 
epoch, and our best fit curve through the data. The size of each point is proportional 
to its weight, i.e., inversely proportional to the square of uncertainty in phase. The error 
bars plotted are ±lcr. From our data through 2005, we obtain a new value for the epoch 
of maximum, T^„^ = 244 2397.9175141 TCB ± 0.41s, a new value for the period, P = 
215.1973888 ± 0.0000004 s, and most importantly, a rate of period change of: 

P = (4.27 ±0.80) X 10-^^s/s. 



^Fitting the whole light curve with a term proportional to sin 



by non-linear least squares 



gives unreliable uncertainty estimates and the alias space in P and P is extremely dense due to 31 yr data 
set span (O'Donoghue 1994; Costa et al. 1999). Our non-linear least squares result is P = (3.38 ± 0.0013) x 
lO-i^s/s. 



- 5 - 



We use linear least squares to make our fit, with each point weighted inversely pro- 
portional to the uncertainty in the time of maxima for each individual run squared. We 
quadratically add an additional Is of uncertainty to the time of maxima for each night to 
account for external uncertainty caused perhaps by the beating of possible small amplitude 
pulsations (Kepler et al. 1995) or the small modulation seen in Figure 1. The amplitude, 
Is, is chosen from the Fourier transform of the (0-C) shown in Figure 2, and is in agree- 
ment with Kopeikin & Potapov (2004) conclusion that the Fourier analysis of the Times 
of Arrival (TOA) gives unbiased information about the noise. Such external uncertainty is 
also consistent with Splaver et al. (2005) who show that the true uncertainties of the times 
of arrival of the milli-second pulsars are generally larger than the formal uncertainties, and 
that a quadratic term is added to them to fit the observations. 



5. Discussion 

While it is true that the period change timescale can be proportional to the cooling 
timescale, it is also possible that other phenomena with shorter timescales can affect P. The 
coohng timescale is the longest possible one. 

As a corollary, if the observed P is low enough to be consistent with evolution, then 
other processes, such as perhaps a magnetic field or diffusion induced changes in the boundary 
layers, are not present at a level sufficient to affect P. 

For the first time we also report on the search for the rates of period changes for the 
other relatively large amplitude modes of G 117-B15A, dP/dt = (36.0 ± 7.2) x 10"^^ s/s 
for the 270s periodicity, and dP/dt = (74 ± 15) x 10"^^ s/s for the 304s periodicity. They 
are much larger than the one derived for the main pulsation, which has a region of period 

formation deeper in the core, but are in line with the measurements for the 274s modes in ZZ 
Ccti (Mukadam et al. 2003). The models to date, without rotation, differential rotation, and 
magnetic fields, do not explain such values, or the chromatic amplitude changes reported by 
Kotak et al. (2004). Kepler et al. (1995) show the 270s and 304s periodicities have significant 
amplitude changes, indicative of multiple components or resonances. 

5.1. Theoretical Estimates and Corrections 

5.1.1. Proper Motion 
Pajdosz (1995) discusses the influence of the proper motion of the star on the measured 
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P: 

-Pobs = -Pevol {l + Vr/c)+ PVr/ C 

where Vr is the radial velocity of the star. Assuming Vj./ c ^ 1 he derived 

Pp^ = 2.430 X 10-^«P[s] {iJi[' /yr]fd[vc] 

where Ppm is the effect of the proper motion on the rate of period change, P is the pulsation 
period, // is the proper motion and d is the distance. The proper motion, // = 0.136 ± 
0.002 "/yr, and the parallax, H = (0.0105 ± 0.004) ", are given by van Altena et al. (1995). 
But the parallax has a large uncertainty and does not agree with other estimates of the 
distance. Munn et al. (2004) measured the same proper motion for G 117-B15A and G 117- 
B15B: = (— 146±2.6) mas/yr, and = (— 1.0±2.6) mas/yr, including 5 USNO-B epochs 
and SDSS positions. Considering = 15.50 ±0.02 (Silvestri et al. 2002), if we use Bergeron 
et al. (2004) estimate of My = 11.70 ± 0.06^ derived from their optical spectra, we obtain a 
distance of 58 ±2 pc, equivalent to a spectroscopic parallax of 0.0172 ±0.0005". But G 117- 
B15A has International Ultraviolet Explorer (lUE) and Hubble Space Telescope (HST) Faint 
Object Spectrograph (FOS) flux calibrated spectra (Koester, AUard, & Vauclair 1994) that 
can also be used to estimate the distance, if we use the evolutionary models of Wood (1995) 
or Panel et al. (2000) to estimate the radius from Tcs and log (7 and fit the observed fluxes 
to a model atmosphere calculated by Detlev Koester, similar to that described in Finley 
et al. (1997), but with an ML2/a=0.6 convection description consistent with Bergeron et 
al. (1995) and Koester & Holberg (2001). The lUE spectra re-calibrated according to the 
New Spectroscopic Image Procession System (NEWSIPS) data reduction by NASA was 
published by Holberg, Barstow, & Burleigh (2003). From the lUE spectra we get a distance 
of d = 59 ± 6 pc. The HST spectra, shown in Figure 3, fits a distance of d = 67 ± 5 pc. 
Table 1 hsts all the available distances to G 117-B15A. 

Taking the average value of proper motion and distance, we estimate 

Pp^ = (7.0 ±2.0) X 10"^S/s 

and 

P = Pobserved " Ppm = (3.57 ± 0.82) X 10"^^ s/s 



^We estimate the uncertainty in the absolute magnitude from the 300 K external uncertainty in Bergeron's 
temperatures. 
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Table 1: G 117-B15A distance determinations 

Method distance (pc) 

Parallax 95 ± 37 

Spectroscopic parallax 58 ± 2 

lUE flux 59 ± 5 

HST flux 67 ± 4 

Seismology from Bradley (1998) 61 

Seismology from Benvenuto et al. (2002) 63 

Mean 67 ± 14 



5.2. Pulsation Models 

We compare the measured value of P with the range of theoretical values derived from 
realistic evolutionary models with C/0 cores subject to ^^mode pulsations in the temper- 
ature range of G 117-B15A. The adiabatic pulsation calculations of Bradley (1996), and 
Brassard et al. (1992, 1993), which allow for mode trapping, give P~ 2 — 7 x 10~-'^^s/s for 
the i — 1, low k oscillation observed. Benvenuto et al. (2004) estimated the theoretical P 
for the three modes of G 117-B15A, even allowing for a dynamical change on the gravity 
constant, as we show in section 1. The observed P/P = 1.9 x 10^ yr, equivalent to 1 s change 
in period in 8.9 million years, is within the theoretical predictions and very close to it. We 
have therefore measured a rate consistent with the evolutionary time scale for this lukewarm 
white dwarf. 



5.2.1. Core Composition 

For a given mass and internal temperature distribution, theoretical models show that 
the rate of period change increases if the mean atomic weight of the core is increased, for 
models which have not yet crystallized in their interiors. As the evolutionary model cools, 
its nucleus crystallizes due to Coulomb interactions between the ions (Lamb & van Horn 
1975), and crystallization slows down the cooling by the release of latent heat. Montgomery 
& Winget (1999) describe the effect of crystaUization on the pulsations of white dwarf stars, 
but G 117-B15A is not cool enough to have a crystaUized core (Winget et al. 1997), or even 
for the convective couphng described by Fontaine et al. (2001) to occur. 

The heavier the particles that compose the nucleus of the white dwarf, the faster it cools. 
The best estimate of mean atomic weight A of the core comes from the comparison of the 
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observed P with values from an evolutionary sequence of white dwarf models. Brassard et al. 
(1992) computed the rates of period changes for 800 evolutionary models with various masses, 
all with carbon cores but differing He/H surface layer masses, obtaining values similar to 
those of Winget et al. (1981), Wood & Winget (1988), and Bradley & Wingct (1991). The 
average value of P for all £ = 1, 2 and 3 modes with periods around 215 s in models with 
an effective temperature around 13,000 K, and a mass of 0.5 Mq, is: P(C core) = (4.3 ± 
0.5) X 10-^5 s/s. Benvenuto et al. (2004) C/0 models give P(C/0 core) = (3-4) x lO-^^s/s. 
Using a Mestel-hke cooling law (Mestel 1952; Kawaler et al. 1986), i.e., T oc A, where A is 
the mean atomic weight in the core, we can write: 

P{A) = (3-4) X 10-^^^s/s. 

The observed rate of period change is therefore consistent with a C or C/0 core. The largest 
uncertainty comes from the models. 



5.2.2. Reflex Motion 

The presence of an orbital companion could contribute to the period change we have 
detected. When a star has an orbital companion, the variation of its line-of-sight position 
with time produces a variation in the time of arrival of the pulsation maxima, by changing 
the light travel time between the star and the observer by reflex motion of the white dwarf 
around the barycenter of the system. Kepler et al. (1991) estimated a contribution to P 
caused by reflex orbital motion of the observed proper motion companion of G 117-B15A in 
their equation (10) as: 

where is the total separation and Mb is the mass of the companion star. In the above 
derivation they have also assumed the orbit to be nearly edge on to give the largest effect 
possible. Only the acceleration term parallel to the line of sight contributes to P. Even 
though G 117-B15A and G 117-B15B are a common proper motion pair (Giclas et al. 
1963; van Altena et al. 2001), there is no other evidence they form a real binary system. 
Kotak et al. (2004) classifies G 117-B15B as an M3Ve from its spectra, obtained with the 
10 m Keck I telescope, and measured loggf ~ 4.5 and Tgff ~ 3400 K. With = 16.1 
(Harrington & Dahn 1980; van Altena et al. 2001) and M^ ~ 10.4 ± 0.9 (Lang 1991) with 
the uncertainty coming from a possible misclassification in the spectral type of G 117-B15B 
from M3V to M4V. Such uncertainty arises from the TiO and CaOH bands seen by Kotak 
et al. (2004) and the (B-V)=1.63 (Harrington & Dahn 1980). Its spectroscopic distance 
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is = 138 ± 47 pc. G 117-B15B is chromospherically active (Kotak et al. 2004) and a 
flare star (Nather & Mukadam 2004). G 117-B15A and B were imaged by SDSS DR3, with 
ua = 15.92, qa = 15.54, = 15.64, za = 15.80, za = 16.06, and ub = 19.38, ge = 
16.95, Tb = 15.46, is = 13.97, Zb = 13.16, which also measure a separation of 13.4". 
They were also observed by 2MASS (Cutri et al. 2003), with Ja = 15.599 ± 0.065, Ha = 
15.614 ±0.12, and Jb = 11.599 ± 0.018, Hb = 11.156 ±0.12, and separation of 13.4" . Using 
Pickles (1998) main sequence spectral templates convolved with SDSS filters, our best fitting 
template is an M4V, with My — 11.54, and a distance of 82 pc. Using Hawley et al. (2002) 
cahbration for the infrared colors, G 117-B15B is consistent with the M3V spectral type, 
and the (i-z) colors correspond to Mf — 10.33, and a distance of 54 parsecs. Using the (i-J) 
colors, Mf = 10.05, we obtain a distance of 61 parsecs. The parallaxes vr"^ = 11 ± 5 mas, 
TT-^ = 5±7 mas (van Altena et al. 2001), and = 4±7 mas (Harrington & Dahn 1980) are 
too uncertain to be used to study the binary nature. Table 2 lists all the available distances 
to G 117-B15B. The mass of an M3.5V should be around (0.30 ± 0.03) M© (Lang 1991). 

Method distance (pc) 

Parallax 200 ± 280 

Spectroscopic parallax 138 ± 47 

SDSS colors 82 

(i-z) colors 54 

(i-J) colors 61 

Mean 107 ± 62 



Table 2: Distance determinations to G 177-B15B 

With a separation around 13.4 arcsec ut — (898 ± 188) AU, corresponding to an orbital 
period around (28 500 ± 2200) years, we get Porbitai < (1-6 ± 1.6) x lO'^^s/s. The large 
uncertainty takes into account the possibility the orbit might be strongly elliptical. Greaves 
(2004), in his section 4.3.1, discusses common proper motion pairs which possibly do not 
form a physical binary. If G 117-B15A and B form a real binary system, the contribution 
of an orbital reflex motion to the observed P might account for one half of the observed P. 

The whole observed rate of period change could also be caused by a planet of Jupiter's 
mass orbiting the white dwarf at a distance of 31 AU, which corresponds to an orbital 
period of 223 yr, or a smaller planet on a closer orbit (see Figure 4). A planet with Jupiter's 
mass any closer to the white dwarf would lead to a larger P. Duncan & Lissauer (1998) 
show that such a planet would survive the post-main sequence mass loss. Note that reflex 
motion produces periodic variations on the O — C, which are distinguishable from parabolic 
variations after a significant portion of the orbit has been covered. 
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As discussed by Damour & Taylor (1991), any relative acceleration of the star with 
respect to the barycenter of the solar system will contribute to the observed P. Their 
equations (2.2) for the differential galactic orbits, decomposed in a planar contribution (2.12), 
where the second term is the proper motion correction, and a perpendicular contribution 
(2.28), applied to G 117-B15A, show the galactic contribution to be exactly the one calculated 
above for proper motion, i.e., the other terms are negligible (2 to 3 orders of magnitude 
smaller) . 

6. Conclusions 

We have measured the rate of change of the main pulsation period for the Tefj ~ 12000 K 
pulsating DA white dwarf G 117-B15A, the first ZZ Ccti to have its evolutionary rate of 
change measured, confirming it is the most stable optical clock known, with a rate of change 
of 1 s in 8.9 million years and a precise laboratory for physics at high energy. It is important 
to notice that mode trapping, a resonance between the local pulsation wavelength with the 
thickness of one of the compositional layers, can reduce the rate of period change by up to a 
factor of 2 (Bradley 1996), but the changes in the trapping layers are still caused by coohng. 

It has taken a huge investment of telescope time to achieve such precision, but not only 
have we measured the cooling rate of this 400 miUion year old white dwarf (Wood 1995), 
excluding the time the star took to reach the white dwarf phase, we have also demonstrated 
it does not harbor planetary bodies similar to Jupiter in mass up to a distance around 30 AU 
from the star or smaller planets with light travel time effects on the white dwarf larger than 
Is. 

We claim that the 215 s periodicity in G 117-B15A is the most stable optical clock 
known. Santra et al. (2004) and Hoyt et al. (2005) discuss projects to build optical atomic 
clocks based on single trapped ions, or several laser cooled neutron atoms, of strontrium or 
ytterbium, expected to reach an accuracy of P < 2 x 10^^^ s/s. Considering their periods 
are 2.5 x 10~^^ s and 1.9 x 10~^^ s, even though they will be more accurate than G 117- 
B15A, they will be much less stable, as their timescales for period changes, P/P, are 125 s 
and 3 h, compared to 2 Gyr for G 117-B15A. Even the Hulse & Taylor's miUisecond pulsar 
(Hulse & Taylor 1975), has a timescale for period change of only 0.35 Gyr (Damour & 
Taylor 1991), but the radio miUisecond pulsar PSR J1713+0747 (Splaver et al. 2005) has 
P = 8.1 X 10-21 s/s, and a timescale of 8 Gyr, and PSR B1885+09=PSR J1857+0943 with 
P = 1.78363 X 10-20 s/s has a stability timescale of 9.5 Gyr (Kaspi et al. 1994). Together 
with these millisecond pulsars, G 117-B15A may be used as the most stable known natural 
frequency and time-keeping standards [e.g. Kopeikin & Potapov (2004)]. The white dwarf 
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has the advantage of not having ghtches and less severe general relativistic corrections. 

It is essential to note that the measured rate of period change is consistent with the 
coohng rate in our white dwarf models only for cores of C or C/0. With the refinement of 
the models, the observed rate of period change can be used to accurately measure the ratio 
of carbon to oxygen in the core of this normal white dwarf, and therefore help in constraining 
the C(q;,7)0 cross section. 

This work was partially supported by grants from CNPq (Brazil), FINEP (Brazil), NSF 
(USA), NASA (USA). This research has made extensive use of NASA's Astrophysical Data 
System Abstract Service. 
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1.5 


2.4 




2443870.946182 


591409 


5.5 


3.1 




2443874.916339 


593003 


2.4 


2.1 




2443959.695117 


627041 


0.1 


2.0 
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Table 3 — Continued 



Time of McLximum 


Epoch of 


(O-C) 


(7 




BJDD 


Mciximum 


(sec) 


(sec) 






698634 


1 6 

X .U 


9 1 

Zi . X 




944'^QQn fifi4fi41 


6'^Q475 


9 7 


1 3 




94441 RQ Q4'SQ^4 


71 1455 


1 


1 6 




Zj'-t'-t'-tZjKj J- . O^^UUU 


7^6998 


-0 7 

U . 1 


9 9 

Z. C/ 




94449S9 SI 8QQ9 


7Q66Q8 

1 OUUC/O 





1 6 

X . U 




944490^^ S'^'^SQfi 


761 1 95 





1 8 

X . O 




9444fiS7 77fi1 74 

ZjH:H:H:UO I.I 1 U X 1 t: 


899985 


5 8 


1 9 

X . C/ 




9444641 694987 


9008*^0 


9 8 


1 1 

X . X 




9444QQ9 7SQ'i'^1 


1041890 


1 


1 6 

X .u 




9444QQ4 fiSQQ'ifi 

^ ^ ^ t? t? . U O t? t? u u 


1 049583 


1 9 


1 9 

X . z 




9444Qqfi 744801 


1 043408 


2 


1 3 

X . (J 




9444QQ7 79 "^640 


1 04*^801 


1 9 


1 9 

X . Z 




944Pin91 71fififi1 


1 05'^4'^4 


1 7 

X . 1 


1 4 




944^170'^ 8finnn4 


1 '^97'^09 


1 9 

X . C/ 


1 7 

X . I 






1 ^^^9668 

X OOC/UUO 


9 4 


1 9 

X . Z 




OAA^T^^ fi4SQ79 


1 340070 

X OtiWW I W 


2 8 


1 3 

X . O 




94461 1 76S716 


1491889 


9 9 


1 9 

X . Z 




9446443 77'i386 


1 694379 


9 8 


1 1 

X . X 




9446468 6301 78 


1 634358 


9 1 


1 3 

X .o 




9446473 71 867Q 


1 636401 


3 


1 6 

X .u 




2446523.620086 


1656436 


2.2 


1.6 




2446524.613917 


1656835 


5.5 


2.5 




2446768.855451 


1754896 


2.9 


1.4 




2446794.935676 


1765367 


2.5 


2.1 




2446796.928219 


1766167 


0.3 


1.6 




2446797.924535 


1766567 


3.1 


1.3 




2446798.903378 


1766960 


2.6 


1.8 




2446823.663537 


1776901 


3.1 


1.9 




2446825.651132 


1777699 


3.7 


1.5 




2447231.328096 


1940575 


3.7 


1.9 
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Table 3 — Continued 



Time of McLximum 


Epoch of 


(O-C) 


(7 




BJDD 


Mciximum 


(sec) 


(sec) 






1 940689 


5 1 

U. X 


3 5 

O.O 




94479'^9 '^Qfifi9fi 


1941004 


5 


1 6 

X .U 




94479'^9 fi9'^9Q1 


1941095 


5 9 


9 9 




94479'?'? '?4'?nQn 


1 Q41 '?84 


4 5 


1 '? 

X . o 




94479'?'? fi'?4^nfi 


1 Q41 501 

X iJ'-t X tj W X 


4 7 


9 '? 
z. o 




94479'?4 '?! 047^; 


1 Q41 77fi 

X C/H: X 1 1 U 


6 8 


'? 9 




94479'?^; '?! '?9Pin 


1 Q491 75 


5 9 


1 4 




94479'?^ fin71 fiS 


1 04990'? 


6 4 


9 1 

Z. X 




94479'?fi fil nQ99 


1 949696 


6 9 


1 6 

X .U 




9447'S8Q '^7'i1 Q8 


9084398 


3 9 


1 4 




94475Q4 '?'?17'?,'i 


2086318 


5 2 


1 6 

X . W 




9447^0"^ 39'^ni8 


9086716 


3 5 


9 




9447PiQfi '?1 1Qn7 


90871 1 '? 

ZjUO I X xo 


10 1 

XU . X 


9 '? 




9447PiQ7 '?1 Pifin9 


908751 6 


4 8 


1 7 

X . 1 




9447PiQS '?! Q'?'?Q 


908701 

ZjUO 1 C/ XC/ 


'? 1 

O. X 


'? 1 

O. X 




24474QQ n72n'?fi 


9048072 


6 5 


3 2 




9447Pi'?9 7fi87QQ 


9061601 

ZjUU XUU X 


1 '? 

X . o 


1 4 




94478'^'^ 84fi'^9'i 


919051 1 

X C/L/U X X 


4 


9 1 

Z . X 




94478'Sfi S'^9fiQ7 


9191710 

J-Cf J- 1 XL/ 


5 9 


1 9 




2447Q1 8 fi44fi'?n 


221 6527 

ZZ X WO i 


2 6 


3 1 

(J. X 




2447920.619811 


2217320 


6.7 


3.3 




2447952.622834 


2230169 


-3.3 


2.9 




2447972.620899 


2238198 


9.6 


6.1 




2447973.709340 


2238635 


9.7 


2.6 




2447973.741682 


2238648 


6.5 


1.4 




2447978.770467 


2240667 


10.0 


2.1 




2447979.781717 


2241073 


11.8 


3.1 




2447980.319627 


2241289 


4.6 


3.5 




2447977.403038 


2240118 


7.5 


2.3 




2447978.327055 


2240489 


4.3 


3.3 
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Table 3 — Continued 



Time of McLximum 
BJDD 


Epoch of 
Mciximum 


(O-C) 

(sec) 


(7 

(sec) 




9447Q7Q "^"SSI 8Q 


9940003 


9 6 


3 4 




9447Q7Q '^'i814'i 


9940Q03 


-1 9 


4 9 




9447078 finiflfiQ 


99405QQ 


7 4 


9 5 




9447Q8n fi9ini7 




5 8 


4 




9447Q8n 7S9Q9Q 


9941475 


7 9 


9 S 




9447Q81 ^^9^01 8 


9941 PiQ">, 


8 4 


1 4 

X . t: 




9447Q81 ^09^0^^ 


9941 800 


5 7 


1 4 

X . t: 




9447081 77Q18^ 


9941875 


4 8 


1 1 

X . X 




9447Q89 39Qfifi3 


99490Q6 


7 4 


1 8 




9447Q89 743nQ3 


99499fi9 


5 


1 9 




9447Q83 7344nn 


2242660 


5 4 


1 2 

X . ^ 




9447Q7Q 9810 "17 


9940879 


9 5 


9 9 




9447080 9948QQ 


99419Pi1 


-9 4 


9 Q 




9447Q84 7'^^fi78 


9943069 


6 5 


1 1 

X . X 




944894^ 794fififi 


9S47847 


O .O 


5 1 




9448967 7QQQS9 


935671 


5 2 


9 3 




9448^^94 fi97Q79 


9*^70596 


4 S 


1 9 




9448395 7n8Q38 


937QQ60 


4 1 


1 3 

X .o 




9448398 5Q39n8 


9381 1 18 


6 4 


1 6 

X .u 




9448331 fifi1735 


9389350 


4 


1 9 

X . Z 




2448238.571479 


2344975 


8.3 


2.2 




2448622.833258 


2499253 


3.3 


1.8 




2448680.642683 


2522463 


6.3 


1.2 




2448687.614155 


2525262 


4.0 


1.2 




2448688.597979 


2525657 


3.4 


1.2 




2449062.660365 


2675840 


4.2 


1.6 




2449063.609354 


2676221 


6.7 


1.9 




2449066.615640 


2677428 


6.5 


1.4 




2449066.371558 


2677330 


7.2 


2.0 




2449066.326737 


2677312 


8.2 


2.6 
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Table 3 — Continued 



Time of McLximum 


Epoch of 


(O-C) 


(7 




BJDD 


Mciximum 


(sec) 


(sec) 






9678593 


fi 4 


1 7 

X . / 




944Q9Q8 9'^Q987 


9770493 


8 5 


4 1 




944Q9Q8 "^04041 




8 9 

O.Z 


4 1 




9440904 '?TA9P,A 


9788807 


5 5 


4 1 




944Q9Q4 9qS8Q7 


97fi88'^0 


-0 5 


4 1 

t:. X 




944Q9QPi 4'^q^S'^ 


97fi0900 


-A 


Pi 1 

U. X 




944Q9QPi 404^^87 


97fiOS91 

Z 1 Uc/OZ X 


O.O 


7 1 

1 . X 




944qn'^fi xnQ9fin 

Z/ ^^t 'rt c/ U O U . O U c/ u u 


9fifi54fi1 


9 4 


9 9 

Z . Z 




944Qn3s fi773nn 


966691 1 


3 1 

O. X 


9 9 




9440040 fiS73in 


9667018 

zuu / uxo 


3 6 


A 1 




9440041 filfiSfin 


9667301 


4 


4 1 




9440700 793S88 


9071765 


5 6 


1 

X .O 




9450497 09n0fin 


3993081 


8 9 

(J . Z 


3 8 




9450490 07S949 


S994805 


9 7 

Z. ( 


9 4 




945n4'^n 01 4770 


S9951 8S 


Fi 


9 5 

Z. tj 






S99555fi 


7 5 


1 5 

X . O 




945n4'^4 01 9*^09 


'?99R788 

OZZU 1 oo 


8 8 

O.O 


9 
z. u 




9450436 090898 


3997508 

OZZ / UC/O 


5 4 


1 7 

X . / 




9450483 6331 80 


'^94fi'^40 


fi 


1 8 

X .o 




9451940 5080060 


3553878 


10 1 

XU. X 


1 3 

X .o 




2451249.7632895 


3553944 


9.7 


1.7 




2451250.6126098 


3554285 


8.7 


2.0 




2451526.8772586 


3665203 


10.4 


1.2 




2451528.8523866 


3665996 


10.0 


1.5 




2451528.9196061 


3666023 


7.4 


1.4 




2451528.9868422 


3666050 


6.3 


1.9 




2451529.8585943 


3666400 


6.6 


2.0 




2451530.9097492 


3666822 


13.4 


2.41 




2451960.8561629 


3839442 


10.1 


1.62 




2451962.7864775 


3840217 


11.3 


1.48 
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Table 3 — Continued 



Time of McLximum 
BJDD 


Epoch of 
Mciximum 


(O-C) 

(sec) 


a 
(sec) 






38491 89 


8 1 


1 93 




94^1 Q88 7Q1Q779 


3850658 


10 5 


9 00 




94^1 QQD 784'i9^'i 


■^851458 


8 8 


1 59 




94'i9n'^7 8479^8'^ 


QS7n97'? 


10 1 








^^87^^489 


19 5 


1 85 






^^045857 


7 


1 ^^4 




94^^999^; Q^QSQ97 


•^045870 


8 

0. U 


fi5 




94'i99fi'^ 88'^481 n 




10 fi 


58 

yj .00 




94 "1931 fi fi4499n'i 


3980791 


1 3 1 


1 




94'S9317 SQQ'i1fi4 


3989988 


19 9 


67 




945931 Q 7QQQ41 7 


3989691 


12 


79 




9459'^1 7 fi47Q75n 




10 


95 




94Pi9S91 S'?4S'^44 




1 1 4 


1 9'? 




94P;9'^99 79fi^9fifi 


3984379 


9 9 


3 03 

xj .\JxJ 




94^9^^! 9 741 9881 


SQ847Sn 


1 1 4 


'\ 5 




94^9^7^^ 8840808 


4n8QQ8'^ 


Q Q 


1 1 

X . X 




9459^7*^ fi8SQ7n9 


4nQn495 


19 fi 


1 9 






4090791 


10 fi 


fi8 




94'S937'i fi3Q97nQ 


41994fi5 


19 8 


1 




941^9374 77nnn7n 

Zi^UZiO / / / yjyjyj 1 \J 


4194076 


19 5 


1 93 




2452581.9494464 


4134940 


12.0 


1.83 




2452583.9095168 


4137288 


11.5 


1.74 




2452584.8812628 


4144503 


13.8 


0.92 




2452585.9821875 


4148953 


12.0 


1.16 




2452586.8937641 


4146104 


12.3 


0.7 




2452665.7845581 


4255806 


13.7 


0.72 




2452669.7970851 


4256260 


10.6 


1.52 




2452696.8561592 


4257768 


12.8 


0.96 




2452724.6624548 


4267394 


10.5 


0.81 




2453381.7442572 


4409917 


11.5 


1.41 
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Fig. 1. — (O-C): (Observed minus Calculated times of maxima) for the 215 s pulsation of 
G 117-B15A. The size of each point is proportional to its weight, i.e., inversely proportional 
to the uncertainty in the time of maxima squared. We show 2a error bars for each point, 
and the line shows our best fit parabola to the data. The error bars plotted are those before 
adding the external uncertainty of Is quadratically, discussed in the text. The fact the line 
does not overlap these error bars arc a demonstration they are underestimate. Note that as 
the period of pulsation is 215.197s, the whole plot shows only ±36deg in phase. 
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Fig. 2. — Fourier transform of the (0-C): (Observed minus Calculated times of maxima) for 
the 215 s pulsation of G 117-B15A, after we subtract the parabola. Kepler's law (1616) shows 
that a* ^ -^73 {^Y^^ P'^^^: considering Mp <C M*, from which we get (light-second) ~ 
0.565^P^/^(yr) i.e., a Is light travel time, as the limit shows, for planets with masses 
between 1.6 Mj and 0.27 Mj=87 M®, for periods between 6 months and 7 yr, respectively. 
The highest peak, at a frequency of 1.85 yr~^, corresponding to a period of 6.5 months, is 
extremely hard to study in our data set, as we normally observe for only 3 months separated 
by 1 year. These variations can be caused by beating of very closely spaced pulsations 
(Kepler et al. (1995)), and we must take into account that multiplets cause non-sinusoidal 
variations, as the triplet in G 226-29 studied by Kepler et al. (1983). 
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Fig. 3. — Fit of a model atmosphere with T^q = 12 000 K, logs' = 8-0, d=67pc, to the HST 
FOS spectra of G 117-B15A. 
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Orbital Separation (AU) 



Fig. 4. — Planet exclusion region around G117-B15A. The filled dots represent known extra- 
solar planets around their stars, while the open circles represent solar system planets around 
the Sun. The short period limit is the limit on the mass of a planet that would produce a 
peak in the Fourier transform smaller than those seen in Fig. 2. The long period limit is the 
mass limit of a companion that would account for the observed P. Any planet above both 
hnes would have been detected. 



